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ABSTRACT 


A theoretical  study  of  the  use  of  a heterodyne  reception  CC^  laser 
radar  for  imaging  and  target  detection  is  reported.  Specifically,  a 
mathematical  system  model  for  the  radar  is  developed,  incorporating  the 
statistical  effects  of  propagation  through  atmospheric  turbulence,  target 
speckle  and  glint,  and  heterodyne -reception  shot  noise.  This  model  is 
used  to  find  the  image  signal-to-noise  ratio  of  a matched-filter  envelope 
detector  receiver  and  the  target-detection  probability  of  the  optimum 
likelihood -ratio  processor.  For  realistic  parameter  values  it  is  shown 
that  turbulence- induced  beam  spreading  and  coherence  loss  may  be  neglected. 
Target  speckle  and  atmospheric  scintillation,  however,  present  serious 
limitations  on  single-frame  imaging  and  target  detection  performance. 
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1 . INTRODUCTION 
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The  proliferation  of  optical  technology  that  has  followed  the  in- 
vention of  the  laser  offers  new  technical  options  for  a variety  of  tactical 
target  detection  and  imaging  scenarios.  In  particular,  both  pulsed  laser 
range-finders  and  imaging  radars  are  promising  candidates  for  immediate 
development  and  deployment  [1],  [2].  The  present  document  constitutes  a 
theoretical  study  of  the  use  of  a heterodyne -reception  CC>2  laser  radar  for 
target  imaging  and  detection.  This  study  had  three  primary  objectives: 

1)  the  development  of  a mathematical  system  model  for  the  CO., 
radar,  incorporating  the  statistical  effects  of  the  propagation 
medium,  the  target  characteristics,  and  the  heterodyne-detection 
process; 

2)  the  use  of  this  system  model  to  investigate  the  influence 
and  interplay  of  propagation  and  target  effects  on  imaging  and 
target  detection  performance; 

3)  the  identification  of  areas  where  there  are  significant 
gaps  in  the  knowledge  needed  to  quantify  the  operational  char- 
acteristics of  the  CO^  radar. 

It  seemed  clear,  a priori,  that  the  entire  area  of  weather  effects, 
i.e.,  the  propagation  characteristics  of  the  atmosphere  under  low-vis- 
ibility conditions,  was  poorly  understood  [3] - [6] . On  the  other  hand, 
the  propagation  effects  encountered  in  clear-weather  conditions,  i.e., 
those  due  to  atmospheric  turbulence,  are  now  well  known  and  adequately 
described  by  relatively  convenient  theoretical  models  [7] -110].  Moreover, 
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reasonable  statistical  models  are  also  available  for  ad  hoc  characterization 
of  optical  targets  [11]- [13].  Thus,  in  addressing  the  preceding  objectives, 
the  system  modeling  and  performance  analysis  was  carried  out  only  for 
clear-weather  propagation  conditions.  As  succeeding  sections  will  testify, 
this  choice  permitted  explicit  derivations  of  the  signal-to-noise  ratio  of 
an  imaging  radar  and  the  receiver  operating  characteristic  (detection 
probability  vs.  false-alarm  probability)  for  target  detection.  Furthermore, 
the  intuitive  understanding  gained  from  the  clear-weather  analysis  has 
proven  valuable  in  preliminary  assessment  of  low-visibility  system  perfor- 
mance. The  remainder  of  this  introductory  section  is  devoted  to  specifying 
the  radar  configuration  that  was  studied,  and  outlining  the  results  that 
were  obtained. 


1.1  CO 2 Laser  Radar  Configuration 

For  the  purposes  of  statistical  system-modeling,  a heterodyne -reception 
CO2  laser  radar  may  be  represented  by  the  block  diagram  of  Fig.  1.  In 
order  to  make  this  block  diagram  conform  with  CO2  systems  under  experimental 
investigation,  the  following  assumptions  will  be  made  at  the  outset. 

The  transmitter  and  receiver  will  be  taken  to  be  co-located  (i.e., 
the  radar  is  monostatic)  with  common  exit/entrance  optics  arbitrarily 
assumed  to  have  an  unobscured  circular  pupil  of  diameter  d in  the  range 
5 cm  ■<  d _<  20  cm.  The  transmitter  laser  will  be  assumed  to  produce  a 
periodic  train  of  rectangular-envelope  purely  sinusoidal  pulses  as  shown 
in  Fig.  2.  The  local  oscillator  will  be  assumed  to  operate  in  a continuous- 
wave  (cw)  mode  producing  an  ideal  monochromatic  wave  displaced  in  frequency 
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TRANSMITTER  ILLUMINATOR 


Fig.  1.  C02  laser  radar  configuration:  (a)  transmitter- to- target  path, 

(b)  target- to- receiver  path. 
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by  an  amount  v^p  (the  intermediate  frequency)  from  the  carrier  frequency 
of  the  transmitter;  the  local  oscillator  power  level  will  be  assumed 
sufficient  for  shot-noise  limited  heterodyne  reception  [14] . 

The  target  of  interest  will  be  located  a distance  L from  the  radar 
along  a line-of-sight  propagation  path  through  the  clear  turbulent  atmos- 
phere; L- values  of  interest  will  lie  in  the  range  1 km  _<  L <_  10  km.  In 
general,  targets  will  be  allowed  to  have  both  specular  ("glint")  and 
diffuse  ("speckle")  reflection  components. 

Mathematically,  the  foregoing  assumptions,  in  concert  with  the  trans- 
mitter pulse  durations  and  pulse  repetition  frequencies  and  target  depths 
forseen  in  realistic  applications  [1],  [2],  have  the  following  implications. 
For  single-pulse  target  detection  and  imaging,  we  may  assume: 

(1)  The  transmitted  wave  has  a monochromatic  linearly  polarized 
electric  field  whose  complex  envelope  is 

Ej,(p)  = (2PT/ceo)1/2  £T(p),  (1.1) 

for  p = (x,y)  a two-dimensional  vector  in  the  transmitter's 
exit  pupil.  In  Eq.  (I . 1)  P^,  is  the  pulse  peak  power,  c is  the 
speed  of  light,  eQ  is  the  vacuum  permittivity,  and  ^ is  the 
normalized  (square- integral  unity  over  the  transmitter  pupil) 
spatial  mode  that  is  transmitted.  Thus,  for  example, 

£p(p)  = (4/7id2)1/2,  (1.2) 

for  | p | j<  d/2  corresponds  to  a collimated  beam  propagating  in 
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the  +z  direction. 


(2)  The  atmospheric  propagation  path  may  be  assumed  to  be  a 
time -independent,  non -depolarizing,  linear  stochastic  spatial 
system  [9],  [10],  [15]- [18].  Thus,  the  complex  envelope  of 
the  electric  field  arriving  from  the  transmitter  in  the  z = L 
plane,  Et(p')  for  p'  = (x',y'),  satisfies 

Ef (p' ) = J dp  Ep(p)  h^Cp'.p).  (1.3) 

In  Eq.  (1.3),  h^  is  the  stochastic  atmospheric  Green's  function. 

(3)  The  target  interaction  may  be  assumed  to  occur  on  the 
plane  z = L so  that  the  complex  envelope  of  the  electric  field 
in  the  z = L plane  that  is  directed  back  towards  the  trans- 
mitter satisfies 


Er(p’)  = Et(p’)  T(p') 


(!.*' 


where  T(p’)  is  the  complex-field  reflection  coefficient  of 
the  target  at  the  point  p'.  This  reflection  coefficient  will 
in  general  be  stochastic  as  it  includes  both  glint  and  speckle 
components . 

(4)  Because  of  atmospheric  reciprocity  [15] , the  complex 
envelope  of  the  electric  field  in  the  receiver  entrance  pupil, 
E^(p) , satisfies 


Ep(p) 


f 


dp'  E (p' ) h. (p'.p). 


(I.: 


Note  that  the  same  Green's  function  appears  in  both  Eqs.  (1.3) 
and  (1.5). 

(5)  Because  of  the  antenna  theorem  for  heterodyne  reception 
[19],  the  local -oscillator  field  and  the  photodetection  process 
may  be  described  as  though  they  were  present  in  the  receiver's 
entrance  pupil.  Thus,  assuming  a local  oscillator  complex 
envelope  E?(p)  given  by 

E^(p)  = (2P^/ce0)1/2  £^(p)  exp(j2TivIFt)  (1.6) 

where  P^  is  the  local -oscillator  power,  and  £^(p)  is  a nor- 
malized spatial  mode,  the  passband- filtered  amplified  photo- 
current (IF  signal)  takes  the  form 

r(t)  = Re  [r(t)exp  (-j2irvIFt)]  (1.7) 

where  the  complex  envelope  r(t)  satisfies  [14] 

r(t)  = (ce0/2 )1/2f  dp  E^p^tp^ntt).  (1.8) 

In  Eq.  (1.8),  n(t)  is  a zero-mean  circulo-complex  Gaussian 
noise  process  whose  power  spectral  density  is  (assuming  a 
rectangular  passband  filter  at  the  IF  frequency) 

Snn(f)  = hvo/n’  for  lfl  - W a-9) 

where  vQ  is  the  transmitter  optical  frequency,  p is  the  quantum 
efficiency  of  the  photodetector,  and  2W  is  the  unilateral  IF 
bandwidth.  Because  the  transmitter  is  actually  pulsed,  Eq.  (1.8) 
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assumes  W > 1/t  . Note  that  for  maximum  heterodyne-detection 

~ _P  . r . . |2  !/2 

efficiency  j^Cp)  = C / dplE^Cp)!  ) is  required.  In  practice, 

however,  £^(p)  = £q*(p)  is  generally  employed,  and  often 
approximates  the  previous  condition. 

The  preceding  mathematical  characterization  of  the  laser  radar  is 

summarized  in  Fig.  3 and  Table  1.  How  the  IF  signal,  r(tj,  is  processed 

depends  on  whether  the  system  is  to  be  used  for  imaging  or  target  detection. 

In  the  former  instance,  r(t)  is  to  be  used  to  estimate  the  average  target 

reflection  strength;  ordinarily  this  will  involve  matched- filter  envelope 

- 1 P 2 

detection,  i.e.,  determining  |t  / r(t)dt|  , followed  by  computer 

P 

enhancement.  (Note  that  a single  photodetector  has  been  assumed  so  that 

a target  image  must  be  built  up  by  scanning  the  transmitter  beam  using  a 

series  of  pulses.  The  extension  of  the  model  of  Eqs.  (1.8),  (1.9)  to  the 

case  of  array  detection  can  be  accomplished,  but  will  be  omitted.)  In  the 

target-detection  application,  r(t)  is  to  be  used  to  determine  the  presence 

or  absence  of  a reflector  in  the  z = L plane;  ordinarily  this  would  be 

P 1 2 

£(t)dt|  with  a threshold  value. 

1.2  Outline  of  the  Report 

The  remainder  of  this  report  is  organized  as  follows.  In  Section  II 
there  is  a tutorial  presentation  of  the  extended  Huygens- Fresnel  principle 
description  of  optical  wave  propagation  through  atmospheric  turbulence. 

Thus  material,  which  comprises  a statistical  characterization  of  the 
Green's  function  h^  appearing  in  Eqs.  (1.3),  (1.5),  will  be  couched  in 
terms  which  emphasize  the  physical  effects  of  turbulence  on  the  fields 


accomplished  by  comparing  |t 
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TRANSMITTER  TARGET-PLANE 

FIELD  FIELD 


Fig.  3.  Suninary  of  CO2  laser  radar  model:  (a)  transmitter- to-target  path, 
(b)  target- to- receiver  path. 


TABLE  I 


SUMMARY  OF  C02  LASER  RADAR  MODEL 
TRANSMITTER  FIELD  Ep(p)  = (2  Rj/  ceQ)1/2  ^,(p)  in  z = 0 plane 

PT  = transmitter  power,  F^Cp)  = normalized  spatial  mode  of  transmitter 

TARGET- PLANE  FIELD  E^p')  - J dp  ET(p)  h^p' , p)  in  z = L plane 

i 

h^Cp',  p)  = atmospheric  Green's  function 
REFLECTED  FIELD  Er(p')  = E^(p')  T(p')  in  z = L plane 
T(p')  = target  complex-field  reflection  coefficient 
RECEIVED  FIELD  E^(p)  ~ f ^p'  Iirtp')  h^Cp'  > p)  in  z = 0 plane 

IF  SIGNAL  r(t)  = Re[_r(t)  exp  (-j2nvjpt)] 

r(t)  = (cCq/2)1^2  J dp  E^p)  F,*(p)  + n(t),  n(t)  = receiver  noise, 


^(p)  = local-oscillator  spatial  mode 
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Et(p’)  and  E^fo).  (These  are  the  fields  which  illuminate  the  target  and  , 

the  receiver's  entrance  pupil,  respectively-)  It  will  be  shown  that  for  I 

the  values  of  d,  L of  interest  the  primary  effect  of  turbulence  on  the  CC^ 
laser  radar  will  be  due  to  scintillation,  i.e.,  beam  spreading  on  E (p') 
will  usually  be  negligible  as  will  be  spatial  coherence  loss  on  E^(p  ) . 

In  Section  III,  a statistical  model  for  the  target  reflection  co- 
efficient T(p')  will  be  discussed.  The  physical  origins  of  the  glint  and 
speckle  components  will  be  described,  as  well  as  their  determination  from 
bidirectional  reflectance  data.  This  determination  is  particularly  import- 
ant, in  view  of  the  wealth  of  bidirectional  reflectance  data  that  has  been 
accumulated  for  various  target  materials  and  shapes  [20]. 

In  Section  IV,  these  results  are  brought  together  to  derive  the  signal- 
to-noise  ratio  (SNR)  obtained  from  matched-filter  envelope  detection. 

A number  of  numerical  examples  are  presented  to  illustrate  the  various 
contributions  of  glint,  speckle,  scintillation,  and  shot  noise  to  image 
SNR.  These  examples  will  indicate  the  importance  of  multi-frame  averaging 
to  obtaining  satisfactory  image  quality. 

In  Section  V,  the  optimum  detection  of  target  presence  or  absence  is 
analyzed.  Results  will  be  presented  for  the  receiver  operating  charact- 
eristic (ROC)  in  the  presence  of  a pure  glint  or  pure  speckle  target.  The 
need  for  further  numerical  work  will  be  indicated. 

The  material  that  will  be  presented  in  Sections  II -V  is  a substantially 
complete  analysis  of  clear-weather  laser  radar  operation.  A suninary  of  the 
key  results  of  this  work  will  be  given  in  Section  VI. 
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II.  CLEAR- ATMOSPHERE  PROPAGATION  MODEL 

■ 

Under  clear-weather  conditions  at  an  operating  wavelength  for  which 

I 

molecular  absorption  is  insignificant,  the  atmosphere  as  an  optical  propa- 
gation medium  differs  from  free  space  in  that  the  former  exhibits  random 
spatio-temporal  refractive  index  fluctuations.  We  shall  refer  to  these 
refractive  index  fluctuations,  which  are  due  to  turbulent  mixing  of  air 
parcels  with  approximately  1°K  temperature  differences,  as  atmospheric 
turbulence.  Although  the  refractive  index  variations  encountered  in  tur- 
bulent air  are  only  a few  parts  in  10^,  their  effect  on  optical  wave  propa- 
gation is  profound. 

At  a qualitative  level,  we  may  gain  an  appreciation  for  turbulence- 
induced  propagation  effects  as  follows.  Let  us  regard  an  instantaneous 
array  of  atmospheric  refractive-index  eddies  as  an  array  of  randomly 
shaped  blobs  with  nominal  spatial  scales  ranging  from  an  inner  scale 
&Q  ^ 10  3 m to  an  outer  scale  LQ  £ 10  - 102  m,  as  shown  in  Fig.  4,  and 
various  refractive -index  values.  Consider  the  effects  of  placing  an  h 
meter  layer  of  such  turbulence  in  the  propagation  path  between  the  radar  > 

r 

and  the  target  in  the  geometry  of  Fig.  3 with  h « L.  Phenomenologically, 
we  have  the  following  situation: 

Transmitter  beam  spread  If  the  turbulent  layer  is  near  the  radar  it  will  1 

randomly  dephase  the  transmitted  field  Ep(p).  In  particular,  when  the 
transmitter  exit  pupil  diameter,  d,  exceeds  the  transverse  phase- coherence 
scale  of  the  turbulence,  the  effective  transmitter  beamwidth  will  be 
turbulence  limited  rather  than  diffraction  limited.  When  the  turbulent 
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TURBULENT  EDDY 


Fig.  4.  Schematic  representation  of  an  array  of  turbulent  eddies  with 
different  refractive  index  values;  tQ  is  the  inner  scale,  Lo  is  the 
outer  scale  of  the  Kolmogorov  inertial  subrange. 
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layer  is  located  close  to  the  target,  there  will  be  essentially  no  trans- 
mitter beam  spreading. 

Target-plane  scintillation  The  random  lensing  of  the  transmitter  beam  by 
the  turbulent  layer  leads  to  constructive  and  destructive  interference 
(amplitude  fluctuations  called  scintillation)  in  the  target-plane  field 
Et(p').  If  the  target  is  smaller  than  the  amplitude-coherence  scale  of  the 
turbulence,  scintillation  merely  modulates  the  reflected  radiant  intensity. 
If  the  target  is  larger  than  this  amplitude -coherence  scale,  scintillation 
presents  itself  as  a speckling  of  the  target  illumination,  |E  (p')|  , which 
the  receiver  may  incorrectly  interpret  as  image  information. 

Receiver  coherence  loss  Because  of  atmospheric  reciprocity  and  the  antenna 
theorem  for  heterodyne  reception,  the  phase  coherence  effects  that  were 
manifest  as  transmitter  beam  spreading  are  also  observed  as  heterodyne- 
receiver  sensitivity  loss.  Specifically,  when  the  receiver  entrance  pupil 
diameter,  d,  exceeds  the  turbulence  phase -coherence  length,  optimum  spatial- 
mode matching  cannot  be  achieved  with  £^(p)  = jj*T(p)  due  to  the  random 
phase  fluctuations  of  the  received  field  E^(p) • This  effect  is  pronounced 
when  the  turbulent  layer  is  near  the  radar,  and  essentially  absent  when  the 
turbulent  layer  is  near  the  target. 

II. 1 The  Extended  Huygens -Fresnel  Principle 

During  the  past  decade,  the  propagation  physics  community  has  taken 
enormous  strides  towards  achieving  a comprehensive  theory  for  wave  propa- 
gation in  the  turbulent  atmosphere  [7] -[10],  [21].  In  this  section,  we 
shall  summarize  the  extended  Huygens -Fresnel  principle  (the  linear- system) 
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model  for  atmospheric  propagation  [9],  [10],  [15],  [17],  [18],  [22] -[25], 
wherein  the  statistics  of  a field  received  from  an  arbitrary  extended 
source  are  obtained,  through  a superposition  integral,  from  the  statistics 
of  a spherical  wave  (point  source) . This  approach  to  the  propagation 
problem  will  quantify  the  qualitative  effects  we  have  just  described  in  a 
form  amenable  to  image  SNR  and  target  detection  ROC  calculations. 

In  the  absence  of  turbulence,  free -space  scalar  paraxial  propagation 
theory  [26]  tells  us  that 

Et(p’)  = f dp  Ej, (p)  (jAL)"1  exp[jkL(l  + | p' -p |2/2L2) ] , (II. 1) 
and 

Ifc(p)  = f dp'  Er(p')  (jAL)'1  exp[jkL(l  + | p' -p| 2/2L2) ] , (II. 2) 

where  A = c/vQ,  k = 2tt/A  are  the  wavelength  and  wavenumber  at  the  trans- 
mitter frequency  vQ.  Equations  (II. 1),  (II, 2)  comprise  the  Huygens -Fresnel 
principle  and  the  reciprocity  theorem  of  Helmholtz;  note  that  they  are  of  the 
forms  assumed  in  (1.3),  (1.5).  In  the  presence  of  turbulence,  the  extended 
Huygens -Fresnel  principle  and  atmospheric  reciprocity  tell  us  that  (1.3), 
(1.5)  are  valid  with 

h^p',  P)  = (jAL)"1  exp[jkL(l  + |p'-p|2/2L2)] 

* exp  (x(p',  p)  ♦ j<Kp',  p) ) » (II. 3) 

where  x and  0 are  real  valued  turbulence  induced  log- amplitude  and  phase 
perturbations  respectively.  Physically,  x(p'»  p)  (p* » p))  is  the  log- 


amplitude  (phase)  perturbation  of  the  field  measured  at  the  point  p'  in 
the  z=L  plane  from  a point  source  located  at  the  point  p in  the  z=0  plane. 
(Via  reciprocity,  the  roles  of  transmitter  and  receiver  may  be  interchanged 
in  the  foregoing  interpretation.)  Mathematically,  it  is  usually  valid  to 
assume  that  x and  $ are  jointly  Gaussian  stochastic  processes  with  known 
means  and  covariances.  For  the  calculations  to  be  made  in  the  sequel,  we 
shall  need  only  the  following  results: 

Wave -structure  function  The  two- source  spherical -wave  wave  structure 
function  defined  by 

D(p'»  p)  = <[x(P1'  + P' * Pj  + p)  - xCp-j^’  , P1)]2> 

+ <[<KPl*  + p' , Pl  + p)  - 4>(p1',  p1)]2>,  (II. 4) 

where  angular  brackets  denote  ensemble  average,  satisfies 

fL 

D(p' , p)  = 2.91  k 2 J dzCn2(z)  (|p’z  + p(L-z)|/L)5/3  (II. 5) 

0 

2 

for  Kolmogorov-spectrum  turbulence  with  Cn  (z)  being  the  turbulence 
strength  profile  along  the  path  from  the  radar  (z  = 0)  to  the  target 
(z  = L) . In  terms  of  this  wave  structure  function  we  have  that  the  Green's 
function  mutual  coherence  function  is 

* . 

^(Pl*  + P'»  Pj  + p)  hj/Pi'’  P1)>  = 

(XL)  2 exp[jk(|Pl'  - P1  + p'  ~ p|2  - 1^'  - px  1 2)/2L) 
x exp  (-D Cp ' t p)/2) . (II. 6) 
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Log -amplitude  covariance  function  The  two-source  log-amplitude  covariance 
function,  defined  by 


C (p' , p)  = <tx(P1'  + p',  Px  + p)  -m  ] [xCp1*  , P-^)  -m  ]>,  CII.7) 


where  m is  the  (assumed)  space -independent  mean  value  of  y (p' , p) , is 

X 


given  in  the  weak  perturbation  regime  by 


cxx&,  p)  - 

4ir^k ^ f dz  J"  du  u Cn^(z)Sn(u)J0(du)sin^[u^z(L-z)/2IcL]  (II. 8) 

0 0 

in  terms  of  an  arbitrary  refractive  index  spectrum  S^(u)  with  d = 

| p*  z + p(L  - z) |/L.  Equation  (II. 8)  has  been  evaluated  when  p = 0 for 

2 

Kolmogorov- spectrum  turbulence  and  a uniform  Cn  profile;  the  resulting 

correlation  length  for  the  log  amplitude  is  approximately  the  Fresnel - 
1/2 

zone  length  (XL)  . The  preceding  weak-perturbation  results  are  valid 
2 - - 

when  a = C (0,  0)  < 0.5;  appreciably  larger  values  of  the  log-amplitude 

X XX 

variance  are  not  encountered  (a  phenomenon  known  as  saturation  of 

2 

scintillation) . It  is  almost  always  reasonable  to  assume  m = -a  , 

X X 

a condition  which  amounts  to  conservation  of  average  irradiance.  For 

2 

Kolmogorov  spectrum  turbulence  with  a uniform  Cn  profile 


o2  = 0.124  C 2 k7/6  L11/6. 
X " 


(II. 9) 


( 

I 

l 

I 

II. 2 Application  to  the  Radar  Configuration 

Let  us  apply  the  model  of  Section  II. 1 to  the  radar  configuration  of 

I 

Fig.  3.  Specifically,  let  us  quantify  the  propagation  phenomena  that  we 
have  described  qualitatively: 

Transmitter  beam  spread  The  average  target-plane  irradiance  that  results 

from  transmission  of  Ep(p)  through  L meters  of  turbulent  air  with  turbulence 
2 

strength  profile  Cn  (z),  0 £ z £ L satisfies 

<|Et(p’)|2>  = <|  J dp  Ej.(p)  hjp’,  p)|2>  (II. 10) 

from  (1.3).  From  (II.3)-(II.6)  we  obtain 

<|Et(p’)|2>  = J dp2  J dp2  ET(pj)  E*(p2)(XL)'2 

x exp[jk( |p^ | 2 - |p2|2)/2L  - jkp'  . (p  1'P2)/L  - D(0,p1-p2)/2] 

(II. 11) 

where  D(0,  p)  = (|p|/p0)5/3  for  pQ  i [2.91  k2^*1  dz  Cn2(z)  (l-z/L)5/3]‘3/5. 

The  quantity  pQ  is  the  turbulence  field  (or  phase)  coherence  length  in  the  * 

transmitter  exit  pupil  plane.  Thus,  when  d « pQ  we  expect  (and,  indeed, 

2 

from  (II. 11)  we  find)  that  <|Et.(p,)j  > is  given  by  the  free-space  target- 
plane  irradiance  pattern.  Furthermore,  when  d > pQ  and  £^(p)  is,  for 
example,  a collimated-beam  mode,  <|Et(p')|  > will  exhibit  severe  turbulence- 
induced  beam  spreading.  Note  that,  in  accord  with  our  intuitive  discussion, 

PQ  weights  turbulence  located  near  the  radar  transmitter  much  more  heavily 
than  it  does  turbulence  located  near  the  target.  It  remains  only  for  us  to 


t 

j 


I 

evaluate  p in  order  to  assess  the  likelihood  of  turbulence -induced  trans- 
o 

mitter  beam  spreading.  In  Fig.  5 we  have  plotted  pQ  vs.  L for  several 
2 

values  of  Cn  assuming  10. 6p  radiation  and  a uniform  turbulence  strength 

2 

profile  along  the  path.  Because  Cn  falls  off  inversely  with  altitude  in 
a fairly  rapid  manner  (below  the  tropopause) , it  would  appear  that  the  CO ^ 
radar  with  5 cm  £ d _<  20  cm  will  only  suffer  from  transmitter  beam  spread- 
ing, over  the  path  lengths  of  interest,  under  conditions  of  strong 
turbulence . 

Target -plane  scintillation  We  shall  assume  henceforth  that  d < pQ  so  that 
for  | p | < d/2  we  may  employ  x(p',  p)  % x(p' » 0),  <Hp\  p)  ^ <t>(p',  0).  We 
thus  find  that  the  target -plane  field  (prior  to  reflection)  is 

Et(p')  ^ (2PT/ceQ)1/2  it(p’)  exp  (x(p',  0)  + j<)>(p',  0)),  (11.12) 

where 

^(p')  e y*  dp  £j.(p)  (j XL) - 1 exp[jkL(l  + |p’  - p|2/2L2)]  (11.13) 

is  the  normalized  free-space  target  plane  field  mode  generated  by  trans- 
mission of  £j.(p).  Equation  (11.12)  shows  the  target-plane  scintillation 
effect  through  the  presence  of  the  x(p’»  0)  term.  For  a uniform  turbulence - 
strength  profile,  10. 6p  wavelength  radiation  will  not  encounter  saturation 
of  scintillation  at  path  lengths  of  interest  except  under  conditions  of 

strong  turbulence  (see  Fig.  6).  Thus,  we  may  generally  employ  the  weak- 

1/2 

perturbation  regime  (XL)  coherence  length  for  x(p'»  0). 
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Fig.  5.  Turbulence  field  coherence  length  p0  vs  propagation  path  length 
L for  conditions  of  weak  turbulence  (C  ^ = 5 x 101'  m L ) moderate 
turbulence  (C^  = 10  ^ m and  strong  (Cn^  = 5 x 10  ^ m tur- 

bulence; 10.6  pm  wavelength  has  been  assumed  throughout. 
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2 

Fig.  6.  Log -amplitude  various  a vs  propagation  path  length  L for  weak 
(Cn^  = 5 x 10  16  m , moderate  (C^  = 10  14  m , and  strong  (Cn^  = 

5 x 10  turbulence;  10.6  pm  wavelength  is  assumed  and  the  weak 

perturbation  theory  is  employed. 


I 


I 

\ 

Receiver  coherence  loss  Assuming  d <<  p„,  there  will  be  negligible  receiver 

o , 

coherence  loss  to  affect  the  heterodyne  reception  of  E^(p).  We  have,  using 
x(p'»  P)  ^ x(p'»  0)  and  <J>(p'»  p)  ^ <Kp' » 0)  f°r  I P I £ d/2,  that 

E^P)  % J dp  Er (p* ) (jXL)'1  exp  [jkL(l  + |p’  - p|2/2L2)] 

x exp  (x(p*  , 6)  + j<j>(p’.  6))  (H.14) 

so  that  the  essential  behavior  of  the  phasefront  of  E^Cp)  is  given  by 

_ 7 

the  spherical  curvature  factor  exp(jk|p|  /2L)  as  would  be  the  case  were  - 

there  no  turbulence  in  the  propagation  path.  Note  that  the  x(p’»  0)  term 
in  (11.14)  represents  additional  scintillation,  beyond  that  included  in  < 

Er(p’)  = Et(p’)  T(p’),  that  is  incurred  on  the  target-to-receiver  propa- 
gation path.  For  the  monostatic  geometry  we  have  assumed,  this  scintilla- 
tion term  is  perfectly  correlated  with  the  scintillation  term  encountered  \ 

in  (11.12). 

III.  TARGET  INTERACTION  MODEL 

The  topics  of  electromagnetic  scattering  from  rough  surfaces,  radar  t J 

target -signature  analysis,  and  optical  speckle  have  a vast  and  growing 
literature  (see,  for  example,  [12],  [20],  [27] - [31 ] ) . Our  purpose  in  the 
present  section  is  to  propose  a statistical  model  for  target  interaction  • 

that  is  analytically  tractable  yet  incorporates  the  accepted  target 
phenomenology  of  glint  plus  speckle  reflections.  We  shall  begin  with  a 
brief  description  of  these  reflection  characteristics. 
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III.l  Planar  Reflection  Model  * 

In  scalar  paraxial  optics,  when  an  optical  beam  which  propagates 
nominally  along  the  +z  axis  encounters  a polished  spherical  surface  of 
intensity  reflection  coefficient  r whose  center  of  curvature  lies  on  or 
near  the  z-axis  and  whose  radius  of  curvature  greatly  exceeds  the  beam 
diameter  of  the  impinging  wave  a planar  reflection  model  may  be  used. 

As  shown  in  Fig.  7,  if  we  designate  a constant- z plane  that  is  rear  the 
reflecting  surface  as  the  input/ output  plane,  this  means  we  can  assume 
that  the  reflected  field  Er(p')  is  given  in  terms  of  the  incident  field 
Et(p’)  by  the  expression 

Er(p’)  = Et(p')  F1/2  exp(-jk| p'  - P'c|2/Rc)  (III.l) 

for  | p* | <_  target  radius  where  Rc  is  the  radius  of  curvature,  p'  is  the 
transverse  location  of  the  center  of  curvature  for  the  surface,  and  we 
have  suppressed  a p'  - independent  phase  factor. 

For  a more  general  polished  reflecting  surface,  with  smoothly  varying 
intensity  reflection  coefficient,  local  radius  and  center  of  curvature, 
the  paraxial  theory  leads  to  an  input/output  relation  of  the  form'*' 

Er(p')  =Et(p')  r1/2(p')  exp(-jk|p*  - p'c(p’)|2/Rc(p’))  (HI. 2) 

*It  would  be  somewhat  more  general  to  allow  for  locally-ellipsoidal  surfaces. 

We  shall  not  do  so  as  (III. 2)  is  sufficient  for  the  physical  arguments  needed 
here  and  in  the  single-glint  approximation. 


i 
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TARGET- PLANE 
FIELD 


Fig.  7.  Geometry  of  the  planar- reflection  model  for  a polished  spherical 
reflector;  the  reflected  field  satisfies  Eq. (III. 1)  where,  for  the  surface 
shown,  R.  > 0. 
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We  shall  assume  that  for  all  targets  of  interest  we  can  assign  a constant-z 
plane  to  be  an  input/output  plane  for  the  target  interaction  in  the  sense 
that  the  incident  and  reflected  fields  are  related  by 

Er(p')  = Et(p')  T(p’).  (HI. 3) 

Equation  (I II. 3)  is  the  model  mentioned  in  Section  1.1.  For  polished 
surfaces  with  smoothly  varying  characteristics,  T,  the  target  field  re- 
flection coefficient,  is  given  explicitly  in  (III. 2).  For  such  targets 
T(p')  characterizes  both  the  reflection  strength  and  surface  shape  of  the 
object  through  its  amplitude  and  phase  behavior.  This  is  not  the  case  for 
rough  targets,  or  for  targets  which  violate  the  paraxiality  assumption 
(e.g.,  in  situations  in  which  shadowing  can  occur).  Nevertheless  in  all 
circumstances  we  shall  regard  T(p')  as  the  target,  i.e. , the  more  the 
radar  can  determine  about  T(p')  the  more  it  has  learned  about  the  target. 

The  planar  reflection  model  Eq.  (III. 3)  has  two  limiting  cases.  The 
first  is  the  case  of  specular  reflection.  This  is  the  limit  wherein  T 
reduces  to  the  form  of  (III. 2)  or  even  (III . 1) . In  such  circumstances, 
the  reflected  field  will  have  directionality  properties  comparable  to 
those  of  the  incident  field;  i.e.,  the  target  acts  like  a mirror  (magnifying/ 
demagnifying  or  distorting  depending  on  whether  (III.l)  or  (III. 2)  holds). 
Under  favorable  alignment  conditions  a specular  object  will  direct  a very 
strong  reflection  towards  the  radar  receiver. 

The  second  limit  of  (III. 3)  is  that  of  diffuse  reflection,  such  as 
will  occur  for  objects  whose  surfaces  are  very  rough  on  spatial  scales 


J 


I 


comparable  to  the  wavelength  of  the  incident  field.  In  such  circumstances, 
the  reflected  field  is  in  essence  spatially  incoherent.  For  laser  beam 
illumination,  however,  the  reflected  field  will  have  temporal  coherence, 
so  that  constructive  and  destructive  interference  effects  will  be  observable 
in  the  irradiance  pattern  at  the  radar  receiver.  This  is  the  well  known 
speckle  effect. 

It  turns  out,  of  course,  that  real  targets  share  aspects  of  both  these 
limiting  cases.  A real  target  may  have  several  "glint"  components,  i.e., 
give  rise  to  several  near-specular  reflections,  in  addition  to  a virtually 
non -directional  speckle  component. 

III. 2 Statistical  Model  for  T(p') 

We  propose  the  following  statistical  model  for  the  target  complex- 
field  reflection  coefficient  T(p')  at  a fixed  wavelength  X.  We  shall 
arbitrarily  distinguish  between  glint  and  speckle  reflections  by  writing 

T(p')  = L,  (S’)  ej6  + T(p').  Cl 1 1 • 4 J 

— -g  — b 

In  Eq.  (III. 4),  T (p* ) denotes  the  glint  reflection  coefficient,  and 
© 

TgCp')  denotes  the  speckle  reflection  coefficient;  these  components  of 

T(p')  will  be  characterized  as  follows. 

Following  standard  practice  [11] -[13],  we  shall  ascribe  T (p')  to  the 

s 

microscopic  surface-height  fluctuations  (surface  roughness)  of  the  target 
and  treat  it  as  a statistical  quantity  with  moments 

<Ts(p’)>  - 0,  (III. 5) 


t 


I 
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^(Pi')  ts(p2’)>  = 0, 


(III. 6) 


<^(Pl’)  t|(p2’)>  = * ^s(p')  fi  Cpj_ * 


(III. 7) 


Equations  (III .5) - (III . 7)  quantify  the  statement  that  a purely  diffuse- 
reflecting  target  turns  a spatially-coherent  illumination  beam  into  a 
spatially- incoherent  reflected  beam;  the  6-function  in  (III. 7)  implies 
this  reflected  beam  has  essentially  no  directionality  to  it.  In  imaging 
a diffuse  target,  the  non-random  quantity  T (p')»  which  is  physically  the 
mean-square  speckle -reflect ion  coefficient  at  location  p',  is  the  target 
information  that  is  sought.  It  turns  out,  however,  that  to  calculate  the 
image  SNR  for  the  diffuse  target  we  must  also  characterize  the  fourth 
moment  of  T (p').  This  is  generally  done  by  assuming  said  moment  can  be 
factored  in  the  manner  that  applies  when  is  a circulo-complex  Gaussian 
process  [32].  Because  we  shall  need  a complete  statistical  characterization 
of  Tg  in  order  to  evaluate  the  ROC  for  diffuse -target  detection,  we  shall 
assume  that  Tg  is  a circulo-complex  Gaussian  process  with  moments  (I I I . 5) - 
(III. 7). 

The  glint  reflection  coefficient,  T (p'),  models  the  specular 

o 

reflection  component  that  is  due  to  the  nominal  (smoothly-varying)  target 

shape.  For  polished  targets  with  smoothly -varying  characteristics, 

T (p')  will  be  given  by  (II I. 2).  For  a "semi-rough"  target  which  has  a 
© 

significant  speckle  reflection  coefficient,  but  still  has  some  glint 
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component,  we  may  continue  to  assume  that  T^(p')  is  of  the  form  (III. 2)  by 

interpreting  r(p')  to  be  the  local  glint  intensity  reflection  coefficient 

at  location  p'  [13).  Such  semi-rough  targets  will  have  |T  (p')|  « T (p'), 

~~s 

although,  because  of  directionality  effects,  the  glint  return  may  dominate 
the  speckle  return  at  the  radar  receiver  entrance  pupil  (see  Section  IV) . 
Note  that  T (p ’ ) will  always  be  treated  as  a deterministic  quantity.  The 
phase  angle  0 that  appears  in  (III. 4)  will  be  taken  to  be  a uniformly 
distributed  random  variable  on  [0,  2tt]  that  is  statistically  independent 
of  the  speckle  reflection  coefficient;  this  phase  angle  models  our  un- 
certainty in  the  location  of  the  target  plane  on  the  scale  of  the  radar 
wavelength.  Although,  for  a purely  specular  target,  we  know  that  both  the 

amplitude  and  phase  of  T (p')  contain  target  information,  we  shall  assume 

S 

2 

the  imaging  radar  is  primarily  interested  in  |Tg(p')|  • For  semi-rough 

targets,  i.e.,  targets  with  non-negligible  glint  and  speckle  reflection 

2 

components,  the  radar  will  be  assumed  to  seek  <|T(p')|  >• 

III. 3 Relationship  to  Bidirectional  Reflectance 

Let  us  examine  the  relationship  between  the  preceding  target 
statistical  model  and  bidirectional  reflectance,  the  target -signature 
quantity  that  is  generally  measured  [2],  [20].  For  the  target  geometry  of 
Fig.  8,  the  bidirectional  reflectance  may  be  defined  as 

P'CA;  fp  fr)  = CA2At)"1<|  f dp'  expUZirf^-fp-p’]  T(p')|2>,  (III. 8) 
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where  Aj,  is  the  target's  projected  area.  In  this  figure,  a target  with 
reflection  coefficient  T(p')  at  wavelength  X is  illuminated  by  a plane  wave 
Et(p')  = exp(j2ir  L • p')  at  this  wavelength,  where  f^  = X * sin0^  x 
(cos<j>^,  sin<^)  gives  the  polar  and  azimuthal  angles  of  incidence.  The  bi- 
directional reflectance  therefore  gives  the  ratio  of  the  average  reflected 

radiance  (W/m  SR)  in  the  direction  ?r  = X i sin6r  (cos^,  sin^)  to  the 

2 

incident  irradiance  (W/m  ).  Equations  (III .4) - (III .8)  yield 

p'(X;  f.;  fr)  = V1/  dP'  TsfP^ 

+ (X2AT)"iy’  dp,'  / dp2'  IgtP]^'  )I^(P2')exp[j  2ir(f'i-fr)  • (px  * -p2' ) ] . (HI. 9) 

Note  from  (III. 9)  that  the  speckle  and  glint  reflection  components  combine 
additively  in  bidirectional  reflectance,  and  that  the  speckle  contribution 
has  no  directionality  whatsoever. 

The  glint  contribution  to  p'(X;  f^;  fr)  can  be  written  as 

(XV1/ dPi '/ dp2’  Ig(Pi’)  Tj(p2’)exp[j21r(fi-fr).(p1'-p2')] 

= (XV1/  dP'  Rgg(p')exp[j2u(fi-fr)*  p'] 

= (X2^)'1  Sgg(fr  - f±),  (III.  10) 

where 

Rgg(Pf)  = f dpx'  T g(px'  + p'/2)  - p 72) 
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is  the  spatial  autocorrelation  integral  of  T , and 

S 


Sgg(f)  = J dp'  Rgg(p')  exp(-j2TT  f . p’) 


is  the  spatial  spectrum  of  T (the  Fourier  transform  of  R ).  Evidently, 

& 56 

in  our  model,  the  directionality  properties  of  p'(A;  fp!  fr)  are  directly 
attributable  to  the  spatial  spectrum  of  the  glint  reflection  component. 

IV.  IMAGE  SIGNAL-TO-NOISE  RATIO 

We  shall  now  bring  the  extended  Huygens -Fresnel  principle  atmospheric 
propagation  model,  and  the  planar  target-reflection  model  to  bear  on 
evaluating  the  performance  of  a C0?  laser  heterodyne  reception  imaging  radar. 

Consider  the  use  of  the  radar  configuration  described  in  Section  1.1 
to  image  a target  at  range  L.  Specifically,  let  us  assume  that  a single 
pulse  is  transmitted  with  spatial  mode 


£T(p)  = (4/iTd2)1^2  exp  Cj 2tt£t* p) , for  |p|  _<  d/2 


(IV. 1) 


i.e.,  the  transmitter  radiates  a collimated  beam  which  propagates  in  the 
direction  of  the  unit  vector 


Lp  = .[27rfT  + (k2  - |2rfT|2)1/2  Iz]/k 


where  Iz  is  a +z  directed  unit  vector.  For  d /XL  < 1,  d < p , as  will 
almost  always  be  the  case  for  target  ranges  of  interest,  far-field  propa- 
gation prevails  with  no  turbulence -induced  beam  spreading.  The  target- 
plane  field,  Et(p'),  therefore  satisfies  (11.12)  with 
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» 

^(p')  * (7tcL2/4)1/2  (j XL)  1 exp ( j kL  + jk|p*  |2/2L) 

xJ1(Trd|p’  - XL£t | / XL) / (wd | p ’ - XLfT(/2XL).  (IV. 2) 

Thus,  the  radar  illuminates  a nominal  XL/d  diameter  region  of  the  target 
centered  on  the  point  p'  = XLfj.  We  shall  assume  perfect  transmitter 
local  oscillator  mode-matching  has  been  achieved,  i.e.,  £^(p)  = £|(p), 
and,  because  d < pQ  has  already  been  assumed,  we  have  that  the  receiver’s 
IF  envelope,  r(t),  satisfies  (1.8)  with 

(ceQ/2 )1//2y*  dp  E^p)  i£*(p)  ^ 

PT1/2/  dp’  It2(p')  T(p’)  exp  (2x(p',  6)  + 2j<|.(p’,  0)).  (IV. 3) 

Equation  (TV. 3)  has  the  merit  of  representing  the  signal  component  of  the 
IF  envelope  r(t)  directly  in  terms  of  a target-plane  integral.  This 
integral  contains  the  interaction  between  the  free-space  radar  beam 

F 

pattern  £t(p'),  the  target's  reflection  coefficient  T(p'),  and  the 
turbulence -induced  propagation  perturbations  exp  (2x(p',  6)  + 2j <j> (p* , 0) ) . 

' t 

IV. 1 SNR  Analysis  i 

Image  SNR  analysis  can  now  proceed  rapidly.  The  radar  receiver's  IF 
envelope  r(t)  consists  of  the  signal  term  (IV. 3),  plus  a statistically 
independent  zero-mean  circulo-complex  Gaussian  noise  process  with  power 
spectral  density  (1.9).  It  is  essentially  straightforward  to  compute  the 
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} 

mean  and  variance  of  the  matched-filter  envelope  detector  output.  For  * 

simplicity,  we  shall  take  W = l/tp,  so  that  the  IF  filter  itself  serves 
as  the  matched  filter  and 

|tp-1y  pr(t)  dt I2  % |r(0)|2. 

0 

With  this  simplification  we  find  that 

<|r(0)|2>  = <|yj2>  + 2hvQ  W/n  (IV. 4) 

and 

Var  (|r(0)|2)=  Var  (|yj2)  + 4 <|yj2>  hv  W/n 

+ (2hv0  W/n)2,  (IV. 5) 

are  the  mean  value  and  variance  of  the  matched-filter  envelope  detector 
output  where 

I = P Tl/Z  f dp’  Ct2(p’)  T(p’)  exp(2x(p’,  0)  + 2j(J>(p’  ,0)) , 

f 

and  expectations  have  been  taken  over  the  turbulence,  target,  and  receiver- 

i -t 

noise  ensembles.  Because  the  2hvQ  W/n  term  in  (IV. 4)  corresponds  to  a 
constant  target- independent  level  (which  is  due  physically  to  receiver 
noise)  we  define  the  image  SNR  to  be 


From  (IV. 4),  (IV. 5)  we  find 


SNR 


(<ll|2>)2 

Var  dill2)  + 4 < | y | 2>  hvQ  W/n  + (2hvQ  W/n)2 


t 

i 


CNR/2 

1 + CNR  Var  (|y;|2)/2  <|;jr|2>2  + (2CNR) 


(IV  .7) 


where  the  carrier-to-noise  ratio  (CNR) 

CNR  = <\y_\2>  / (2hvo  W/n)  (IV. 8) 

2 

is  the  ratio  of  the  mean  signal  and  mean  noise  contributions  to  jr(0) | . 

The  three  terms  in  the  denominator  of  (IV. 7)  represent,  physically, 
the  fluctuations  in  |r (0) ( that  are  due,  respectively,  to  signal  x noise 
fluctuations,  signal  x signal  fluctuations,  and  noise  x noise  fluctuations.  In 
the  absence  of  turbulence  and  target  fluctuations,  we  have  the  Var  (|^J  ) = 

0 (see  below);  at  high  CNR  values  we  then  find  SNR  'v  CNR  / 2.  For 
CNR  _>  5 the  noise  x noise  term  in  the  denominator  of  (IV. 7)  becomes 
insignificant  and  we  can  use 


SNR  £ 


CNR/ 2 

1 + CNR  Var  (|x|2)/2  <lz|2>2 


which  is  equivalent  to  the  universal  curve 


SNR  £ 
SNRsat 


CNR  / 2 SNRgAT 
1 + CNR  / 2 SNRsat 


(IV. 9) 


(IV. 10) 


i 
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in  terms  of  a saturation  signal -to -noise  ratio 


SNRsat  = < IxJ 2>2/Var  (|Z|2).  (IV. 11) 

Equation  (IV. 10)  has  been  plotted  in  Fig.  9;  evidently  increasing  CNR 

beyond  2SNRgA^.  does  little  to  improve  image  signal -to -noise  ratio. 

The  preceding  discussion  implies  that  we  should  focus  our  attention 

on  evaluating  CNR  and  SNR^^,  for  various  target  and  turbulence  conditions. 

2 2 

To  do  this  we  must  find  <|^|  > and  Var  (|yj  ). 

Using  the  results  of  Sections  II  and  III  we  find  that 

<W2>  = PT  A2  e4ox2/  dp'  |£t2(p’)i2  Ts(p') 

2 

* PT  e4°X  / <V  / <ip2'  IJCP21) 


x exp  (-2  Dfpj*  - P2',  0)) 

= PT  e4ax2  [X2/  dp'  |it2(p’)|2  Ts(p*) 

+ / d'P'  Rgg(P,;)  exP[‘2(|P’l  / P0')5/3]]  (IV. 12) 


where,  in  the  second  equality,  R'  is  the  spatial  autocorrelation  integral 
of  it2  Ig  (cf-  (HI-10))  and  po’  = [2.91  k 2 f dz  Cn2(z)  (z/L)5/3]  is 
the  turbulence  field  coherence  length  in  the  target  plane  (cf.  (II. 11)). 
Equation  (IV.  12)  applies  to  an  arbitrary  target  (polished,  rough,  or 
semi-rough)  and  arbitrary  turbulence  levels  (up  to  the  point  d = pQ). 
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Weak  perturbation  turbulence  theory  can  be  used,  in  conjunction  with  our 

target  model,  to  obtain  a general  expression  for  Var  (|yj  ).  We  shall 

not  produce  such  an  expression;  it  is  near-hopeless  in  its  complication. 

2 

Instead,  we  shall  develop  from  (IV. 12)  simpler  results  for  <|^J  > and 
2 

Var  (|y|  ) when  the  target  has  but  a single  glint  (in  addition  to  its 
speckle  component) . 

IV. 2 The  Single-Glint  Approximation 

Consider  the  spatial  autocorrelation  integral 


Rgg(^3  =f  ‘V  + ~p,/2)  + p'n) 


x Ei2&i  - P'/2)  T*^'  - p*/2)  | 


(IV. 13) 


which  enters  into  (IV. 12).  Because  |Ct(p')|2  ^ 0 for  |p'  - ALfT|  > AL/2d 
(recall  (IV.2)),  |R'  (p')|  ^ 0 for  |p'|  > AL/d.  Thus,  if  AL/d  < p '/2 

oo  ^ 


then  the  glint  contribution  to  < |^|  L>  satisfies 

/dp*  ^g(p')  exp[-2(|p*  |/p0')5/3]  ^ / dp’  Rgg(p'), 


(IV. 14) 


i.e.,  the  turbulence  causes  no  beam  spreading  of  the  glint  return.  Un- 
fortunately, AL/d  < pq'/2  is  an  unduly  restrictive  condition,  viz,  it  is 
violated  for  a variety  of  d,  L,  Cn  values  of  interest.  A more  careful 
examination  of  R'  (p')  reveals  that  there  is  a much  weaker  condition  which 

Oo 

permits  us  to  use  (IV. 14).  Indeed,  this  condition,  developed  below,  per- 
mits a useful  simplification  of  the  glint  contribution  to  y_  itself,  before 
ensemble  averaging . 
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Suppose  Tfp' ) is  given  by 

O 

I^(p')  * r1/2  exp(-jk|  p'  - p’J2^  (IV. 15) 

for  | p'  - XLfT|  £ AL/2d,  i.e.,  the  glint  reflection  behaves  like  a spherical 
reflector  over  the  region  illuminated  by  the  radar.  Assuming  R,  « L (which 
should  often  if  not  always  be  the  case),  (IV. 15)  implies  that  |R'  (p')|  % 0 

oo 

dR 

for  ) o' | > c . It  then  follows  that  (IV. 14)  may  be  employed  because 

I* 

d < p , R <<  L have  been  assumed  and  p ^ p ' for  horizontal  or  near-horiz- 
o c o o 

2 

ontal  paths  (uniform  Cn  profiles) . 

We  can  show,  by  direct  consideration  of 

PT1/2/  dp'  it2(p')  Tg(p')  exp(2x(p’,  0)  + 2 j 4> (p * , 0)  + je) 

which  is  the  glint  contribution  to  £,  that  when  (IV. 15)  is  satisfied,  only 

1/2 

a region  of  nominal  diameter  (ARC)  within  the  illuminated  portion  of  the 

target  makes  an  appreciable  contribution.  This  is  the  stationary-phase  region 

[33]  on  the  target;  it  is  the  Fresnel  zone  of  the  curvature  Rc  centered  on 

the  reflection  point  which  produces  a wavefront  directed  towards  the  radar 

receiver  when  illuminated  by  £t(p'),  see  Fig.  10.  For  reasonable  Rc  values 
1/2 

(ARc)  ' <00'  can  be  assumed  hence 

PT 1/2 f dp'  it2(p’)  Tgtp')  exp(2x(p' , 0)  + 2j<J>(p\  0)  + je) 

% PT1/2/  dp*  ^.2(p')  Tg(p')  exp(2x(p’g,  6)  + 2j0(p'g,  0)  + je)  (IV. 16) 
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18-5  -9266 


LOCALLY  - SPHERICAL 
REFLECTING  SURFACE 


z = L 


Fig.  10.  Geometry  for  the  single-glint  approximation;  the  target-plane 
field  has  direction  of  propagation  ip  and  illuminates  a region  of  the 
target  with  nominal  diameter  XL/d  centered  on  the  point  ALfp  in  the  z = L 
plane  (X?p  is  the  projection  of  Ip  on  the  z = L plane) ; over  this  region 
the  target  is  assumed  to  be  a spherical  reflector  with  center  of  curva- 
ture pc’(ALfp)  and  radius  of  curvature  RcfALfp)  (Rq  > 0 is  shown);  the 
only  appreciable  contribution  to  the  received  field  at  z = 0 comes  from 

1/2 

a Fresnel-zone  region  (nominal  radius  (AR^.)  ) about  the  glint  reflec- 

tion point  p ' ; p ' is  found  by  projecting  p ' (ALfj)  onto  the  z = L 

plane  in  the  direction  of  -Lp  for  Rc  < 0 and  Lj,  for  Rc  » 0. 


where  p'  is  the  glint  reflection  point.  (Use  of  (IV. 16)  reproduces  the 
S 

weaker  condition  (IV. 14),  but  the  former  has  additional  advantages  in 
finding  Var(|^|2)). 

We  can  now  state  the  single-glint  approximation.  If,  for  | p * - ALfT| 

£ AL/2d,  T (p')  has  only  a single -stationary  phase  point  for  specular  reflec- 
S 

tion  of  the  beam  mode  j^.(p')  back  towards  the  radar,  and  if  the  Fresnel  zone 

1/2 

(effective  radiating  region)  for  this  point  has  nominal  diameter  (AR^) 

with  R <<  L,  (AR  )^2  < p ' then 
c c o 

z * PT1/2  / dp*  £t2(p’)  T^p')  exp(2x(p',  6)  + 2j«Kp*,  0)) 

+ P T1/2/  dp’  it2(p')  Tg(p')  exp(2x(p'g,  0)  + 2j<j,(p'g,  0)  + je),  (IV.  17) 


where  p'  is  the  glint  reflection  point,  is  the  signal  contribution  to  the 
© 

receiver's  IF  complex  envelope  r (t) . 

From  (IV. 17),  and  the  results  of  Sections  II  and  III,  we  can  show  that 


<iii2> 


% PTe 


4ox 


[A 


lf  & lit 


2(P’)|2 


Ts(p’) 


•/d'p' 


(IV. 18) 


For  an  object  that  lies  entirely  within  the  target-plane  region  illuminated  by 
the  radar,  we  have 

R’g(p')  * lit2(ALfT)|2  exp(j47rfT.p')  Rgg(p'),  (IV. 19) 
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2 

and  the  glint  contribution  to  <|^|  > is  directly  related  to  the  bidirectional 
reflectance.  Using  the  weak-perturbation  turbulence  theory  and  the  single- 
glint  approximation  we  now  get  this  "simple"  expression  for  Var(|yJ  ) 

Var  (|yj2)  % 

2 

A2PT2e8aX  f <V  f dp2'  TgCpj*)  Ts(p2')  |^.2Cp1’)|2 
x!it2(P2’)|2  [2exp(16  C^'  - P2'»  6))  - 1] 

+ 2A2PT2e8ax  f dpx'  Tg (Pj_ * ) llt2(P1,)|2 

x[2exp(16  C - pg\  6))  - 1]  J dp2’  RggCP2') 

♦ PT2e8°x  ( / dp*  RggCp’))2  (e16°X  - 1)  (IV. 20) 

Equations  (IV. 18),  (IV. 20)  will  be  used  henceforth;  in  conjunction  with 
Equations  (IV. 7)-(IV.ll)  they  complete  the  general  SNR  analysis  for  the  imaging 
radar  in  the  single-glint  approximation. 

IV. 3 SNR  Examples 

To  illustrate  the  interplay  between  receiver  noise,  turbulence- induced 
scintillation,  and  target  speckle  in  determining  image  SNR,  we  shall  consider 
a series  of  simple  examples.  At  times,  rather  gross  approximations  will  be 
made  to  arrive  quickly  at  numerical  results.  These  approximations,  however, 
will  not  detract  from  our  main  purpose,  viz,  to  gain  an  intuitive  understand- 
ing of  the  various  SNR  degradations. 
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Case  1 No  turbulence,  specular  target. 


Consider  the  imaging  of  a specular  target,  i.e.,  one  for  which 

T(p’)  = T (p’)  eJ0,  in  the  absence  of  turbulence.  (This  is  the  most  favor- 
S 

able  situation  inasmuch  as  image  signal-to-noise  ratio  will  be  limited  only 
by  receiver  noise.)  We  have  that  the  fundamental  performance  equations 
reduce  to  the  following  set  for  this  case: 

CNR0  /2 

SNR°  = 8 (IV.  21) 

8 1 + (2CNR0  ) 

8 

CNR°g  = (nPT/2hvoW) J dp'  R^g(p’),  (IV. 22) 

where  the  superscript  "o"  denotes  free  space  propagation  and  the  subscript 
"g"  denotes  glint  target. 

For  simplicity,  in  this  example  and  those  that  follow  we  shall  use 
(IV. 19)  for  R’  (p')»  whence  (with  the  aid  of  (III. 10))  (IV. 22)  reduces  to 

oo 

CNR°g  = (nPT/2hvQW)  |£t2(XLfT)|2A2ATp'(X;fT;-fT).  (IV. 23) 

It  is  reassuring  to  see  that  our  theory  reduces  to  such  an  intuitively 
obvious  result.  Equation  (TV. 23)  shows  that  the  no-turbulence  no-speckle 
carrier-to-noise  ratio  (for  a target  smaller  than  the  radar  beam  size) 
equals  the  on-target  irradiance,  |£t(XLf^,) | , times  the  bidirectional 
reflectance  for  specular  reflection  back  towards  the  radar  p' (X;fT;-fT) , 
times  the  target's  area,  Aj,,  times  the  heterodyne  receiver's  antenna  gain, 
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2 9 

A |Ct(ALfT)|  , divided  by  the  locjal -oscillator  shot  noise  power  in  the  IF 
bandwidth,  2hvQW/n. 

Equation  (IV. 23)  is  plotted  in  Figure  11  vs.  path  length  L,  assuming 

the  parameter  values  given  in  Table  2.  Because  the  extended  Huygen-Fresnel 

principle  we  have  employed  neglects  absorption,  (IV. 23)  is  optimistic.  We 

can  correct  for  absorption  in  all  of  our  signal -to-noise,  carrier-to-noise 

results  by  replacing  PT  by  PT  exp (-2  a L)  where  a is  the  10.6pn  absorption 

coefficient;  two  such  absorption-corrected  CNR°  curves  are  included  in 

S 

Fig.  11.  Inasmuch  as  CNR°  _>  5 will  prevail  under  the  conditions  assumed 

o 

in  Fig.  11  we  may  use  SNR°g  % CNRy2  in  lieu  of  (IV. 21).  Note  that  this 
corresponds  to  a curve  of  the  form  (IV. 10),  with  a free-space  propagation/ 
glint  target  saturation  signal-to-noise  ratio  SNR°gSAT  = °°,  as  expected 
from  (IV. 11)  with  Var(|^|2)  = o. 

Case  2 No  turbulence,  diffuse  target 

In  the  absence  of  turbulence  with  T(p')  = T^fp') , we  obtain  the  image 
signal-to-noise  ratio 


SNR°  = 
s 


CNR°  /2 
s 


1 + CNR°s/2  + (2CNR°s)' 


(IV.  24) 


where 


CNR°S  = (nPT/2hvQW)  A 2 f dp'  |it2(p')|2  TsCp ’),  (IV. 25) 

and  the  subscript  ”s"  denotes  speckle  target.  It  follows  from  (IV. 24)  that 

SNR0  increases  monotonically  with  increasing  CNR0  . When  CNR0  > 5, 

5 s s — 
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TABLE  2 


PARAMETER  VALUES  FOR  CNR°  EVALUATION  l ! 

g 

TRANSMITTER  PEAK  POWER  P?  = lO^f  1 

DETECTOR  QUANTUM  EFFICIENCY  n = 0.2 

TRANSMITTER/RECEIVER  PUPIL  DIAMETER  d = 0.1  m j 

7 * 

UNILATERAL  IF  BANDWIDTH  2W  - 2.0  x 10  Hz- 

-20 

CO-  LASER  PHOTON  ENERGY  hv  = 1.87  x 10  u J 
2 o 

GLINT  BIDIRECTIONAL  REFLECTANCE  (see  [2,  Fig.  16])  < 

p* (X;  fT;  -fT)  = 0.S  SR'1 

TARGET  AREA  Ap  = 4(AL/d)  2/tt  I 

ABSORPTION  COEFFICIENTS  (see  [2,  TABLE  2]) 
a - 0 (NO  ABSORPTION) 
a = 0.23  km"1  (1  dB/km) 
a = 0.46  km  1 (2  dB/km) 

■y 
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however,  Eq.  (IV. 24)  can  be  replaced  with 

CNR0  /2 

SNR0  = ^ (IV.  26) 

s 1 + CNR°s  / 2 

which  is  of  the  form  (IV. 10)  with  a free-space  propagation/ speckle  target 

saturation  signal-to-noise  ratio  SNR°s  = 1.  Thus,  in  the  limit 

CNR°  » 1,  we  find  SNR0  £ 1,  which  is  the  well  known  unity  signal-to-noise 
s s 

ratio  of  laser  speckle  111],  [12],  Multiframe  averaging  is  needed  to  overcome 
this  speckle- induced  performance  limit. 

Assuming  here  and  in  subsequent  examples  that  Ts(p')  does  not  vary 
appreciably  over  the  target  region  illuminated  by  the  radar,  (IV. 25)  becomes 

CNR°s  % (qPT/2hv0W)  X2  Ts(ALfT)  f dp'  |it2(p’)|2 

% (nPT/2hvQW)  (d/L)2  Ts(ALfT).  (IV. 27) 

Equation  (IV. 27)  has  been  plotted  vs.  path  length  in  Figure  12,  using 
Tg  = 0.1  and  other  parameters  as  in  Table  2.  (We  have  also  included  in 
Figure  12,  CNR°S  plots  for  two  representative  values  of  absorption  coefficient.) 
For  this  example,  CNR°s  » 1 will  prevail  (except  for  the  case  of  2dB/km 
absorption  at  path  lengths  in  excess  of  6 km),  so  that  image  signal-to-noise 
will  indeed  be  speckle  limited. 

Case  5 No  turbulence,  semi-rough  target 

In  the  absence  of  turbulence  with  a semi-rough  target,  image  signal- 
to-noise  ratio  satisfies 
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SNR°  = 2®72 (IV.  28) 

1 + CNR°s/2  + (2CNR°)  1 

where  CNR°  = CNR°  ; + CNR°  . For  the  parameters  used  in  Figures  11,  12  we 
8 ® 

have  CNR°  >_  5,  so  that  (IV. 28)  can  be  reduced  to 

SNR0  = CNRV2 (IV.  29) 

1 + CNR°/2  (CNR°/CNR°s) 

Equation  (IV. 29)  can  be  put  in  the  standard  form  (IV. 10)  by  identifying  the 
free-space  propagation/ semi -rough  target  saturation  signal-to-noise  ratio 


SNR0  = CNR°/CNR°  = 1 + CNR°  /CNR°  . 
SAI  s g s 


(IV. 30) 


For  the  parameter  values  assumed  in  Table  2,  Figs.  11,  12  tell  us  that 

CNR°g  > CNR°s  > 5 hence  SNR°  'v  SNR°g^y  In  this  case,  image  signal-to-noise 

is  limited  by  the  speckle-image  "noise"  that  appears  superimposed  on  the 

specular  image.  Other  limiting  cases  of  (IV. 29)  are  also  possible,  e.g., 

CNR0  » CNR°  » 1 which  reduces  SNR°  to  SNR°  etc. 
s g s 

Case  4 Turbulence,  specular  target 

In  imaging  a specular  target  through  turbulence  we  find  that  image 

signal-to-noise  ratio  satisfies 


SNR  = 
g 


CNR  / 2 




1 + CNR  (e16ox  - l)/2  + (2CNR j"1 
g g 


(IV. 31) 


* 

where  CNR  = CNR°  e^°x.  Assuming  CNR  5,  (IV. 31)  takes  the  standard 

O g g 

form 
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CNR  /2 

SNR  = 8 

8 1 + CNR  /2SNR  cvr 

g gSAT 


(IV. 32) 


SNRgSAT  ‘ Ce16°X  -W1- 


(IV. 33) 


We  see  that  for  a > 0,  SNR  initially  increases  with  increasing  CNR  , but 
X 8 & 

reaches  the  scintillation-limited  value  (IV.  33)  once  CNR^  >>  2SNRgSAT' 

Evidently,  turbulence -induced  scintillation  on  a specular -target  image  is 

like  the  speckle  fluctuations  of  the  free-space  image  of  a diffuse  target, 

i.e.,  high  carrier-to-noise  ratio  does  not  lead  to  high  image  signal-to-noise 

ratio.  Multiframe  averaging  will  be  required  to  overcome  this  limit. 

To  gain  some  quantitative  appreciation  for  the  foregoing  scintillation- 

2 

limited  signal-to-noise  ratio,  we  have  plotted  SNRgg^,  vs.  a in  Fig.  13,  and 

2 2 
o^  vs.  path  length  L (assuming  a uniform  Cn  profile)  in  Fig.  14.  These 

figures  show:  for  L ^ 1 km  and  realistic  C2  values  (10  ^ -10  ^ m~ 2^) 

2 4 2 2 

a may  range  from  less  than  10  ^ to  in  excess  of  10  , a values  greater  than 

X X 

10  2 imply  SNRgSAT  _<  60. 

Case  5 Turbulence,  diffuse  target 

For  a diffuse  target,  image  signal-to-noise  ratio  in  the  presence  of 
turbulence  satisfies 


SNR  = 
s 


CNR  /2 
s 


1 + CNRs[l  + 2(e16oX  - 1)  C ]/2  + (2CNR,.)'1 


(IV. 34) 


MM 


Fig.  14.  Profiles  of  constant  log-amplitude  variance  ax  vs  path  length 
L and  turbulence  strength  C 2;  10.6  ym  wavelength  is  assumed. 


In  (IV.  34),  CNRg  = CNR°s  ex  and 


/dpi'/dpz’  |jt2(p1,)|2  ljt2(P2,)|2[exp(16  Cyy(p  '-P2,>0)-1] 
(e16ox  - 1)  ( / dp  |^2(p’)|2) 


(IV. 35) 


is  the  log-amplitude  aperture  averaging  factor  [7],  [34].  It  is  easily  shown 
that  5 _<  1.  When  AL/d  < (XL)1/2  so  that  the  radar  illuminates  a single  log- 
amplitude  coherence  area  on  the  target  plane  we  have  C ^ 1;  when  AL/d  » (AL)^2 
we  have  approximately  AL/d  uncorrelated  log-amplitude  regions  illuminated  by 

r\j  2 

the  radar,  whence  c ^ d /AL.  Both  of  these  regimes  can  be  conveniently 
summarized  by  using 


r ^ d2/AL 

C ^ * 

1 + d / AL 


(IV.  36) 


in  place  of  (IV. 35).  (This  result,  remember,  presumes  no  transmitter  beam 
spreading  due  to  turbulence,  and  uses  the  weak-perturbation  log-amplitude 
coherence  length . ) 

When  CNRg  5,  (IV. 34)  reduces  to 

CNR/ 2 

SNR  = 2 (IV.  37) 

1 ♦ CNRs/2SNRsSAT 


where 


SNR 


‘sSAT 


= [1  + 2 (e16ox  - DC]"1 


(IV.  38) 


is  the  atmospheric  propagation/ speckle  target  saturation  signal -to-noise  ratio. 
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In  Fig.  15  we  have  plotted  SNRgc^j  vs.  L for  three  representative  Cr  values. 

These  curves  are  obtained  assuming  a uniform  Cn  distribution,  d = 10  cm,  and 

a target  that  fills  the  main  lobe  of  the  illuminator  beam  in  the  plane  z = L; 

the  curves  are  cut  off  at  path  lengths  before  saturation  of  scintillation  is 

encountered.  Note  from  Figure  15  that  SNRsg^  < 1;  SNR^^-j,  = 1 corresponds, 

physically  to  a signal-to-noise  ratio  limited  by  target  speckle,  the  presence 

2 -12-2/5 

of  strong  turbulence  (C  =10  m ) or  a long  path  length  (L  'v  10  km) 

2 -14-2/5 

with  moderate  turbulence  (Cfi  =10  m ' ) may  significantly  reduce  SNR^^-j, 
from  the  speckle -limited  value  of  1.  To  illustrate  the  aperture  averaging 
effect  in  Figure  15,  we  have  plotted  c vs.  L for  d = 10  cm  in  Figure  16, 
assuming  the  target  fills  the  main  lobe  of  the  illuminator  beam.  Without 
aperture  averaging,  the  SNRsg^j  values  in  Figure  15  would  be  more  than  one 
order  of  magnitude  lower  at  L = 10  km.  In  this  regard,  it  is  worth  noting 
that  targets  which  do  not  fill  the  main  lobe  of  the  illuminator  beam  do  not 
achieve  the  full  aperture  averaging  predicted  by  (IV. 36). 

V.  RECEIVER  OPERATING  CHARACTERISTIC  FOR  TARGET  DETECTION 

The  single-pulse  target-detection  problem  for  the  CO2  laser  heterodyne 
reception  optical  radar  is  as  follows.  There  are  two  possible  hypotheses 
(target  present  or  target  absent  at  transverse  coordinates  p’  = ALf-j,  in 
the  plane  z = L) . If  there  is  no  target  at  p'  = ALfp  in  the  z =L  plane 
(hypothesis  Hq),  then  the  receiver's  IF  complex  envelope,  r(t),  satisfies. 


r(t)  = n(t)  for  0 < t < t 


(V.l) 
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>4  -2/31 


10.6  pm 

0.1  m 


18-5-9531 


Fig.  15.  Atmospheric  propagation/ speckle  target  saturation  signal -to-noise 
ratio  SNRsg^p  vs  path  length  L for  various  turbulence- strength  (C  2)  values; 

10.6  pm  wavelength,  d = 0.1  m,  weak-perturbation  propagation  theory  and  a 
target  area  which  fills  the  main  lobe  of  the  illuminator  beam  have  been 
assumed. 


■ 


10  102  103  TO4 

L ( m ) 


. 

Fig.  16.  Aperture  averaging  factor  ? vs  path  length  L for  a 10  cm  pupil 
diameter  at  10.6  jim  wavelength. 
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where  n(t)  is  the  complex-Gaussian  noise  process  with  power  spectrum  (1.9) 
and  W > t has  been  assumed.  If  there  is  a target  at  p'  = ALfp  in  the 
z = L plane  (hypothesis  H^) , then  the  receiver's  IF  complex  envelope  obeys 


r(t)  = + n(t)  for  0 £ t < t (V.2) 

where  (as  in  Section  IV) 

Lm  P TV2/  dp'  it2(p')  T(p')  exp(2x(p' , 0)  + 2j  <p(p> , 0))  (V.3) 

Our  objective  is  to  process  the  data  (r(t)  : 0 £ t <_  t^}  in  such  a way  that 
the  conditional  probability  we  say  there  is  a target  present  given  there  is 
one  there  is  maximized,  subject  to  the  constraint  that  the  conditional 
probability  we  say  a target  is  present  given  there  is  no  target  there  does 
not  exceed  a specified  level.  This  is  the  well  known  Neyman- Pearson  approach 
to  binary  hypothesis  testing  [35] . The  optimum  decision  rule  (data  processor) 
is  the  likelihood  ratio  test  (LRT)  whose  threshold  level  is  chosen  to  satisfy 
the  false-alarm  probability  (probability  you  say  given  Hq  is  true)  con- 
straint with  equality.  The  performance  attainable  with  an  LRT  is  conveniently 
summarized  by  its  receiver  operating  characteristic  (ROC) , which  is  a para- 
metric plot  of  the  detection  probability  (probability  you  say  given 
is  true)  vs.  the  false-alarm  probability  Pp  obtained  by  varying  the  threshold 
value. 

V.l  The  Likelihood-Ratio  Test 

The  well  known  results  for  signal  detection  in  additive  white  Gaussian 


f 1 

« 


noise  [36]  permit  us  to  express  the  LRT  for  the  C02  radar  in  the  following 


form  when  W = 1/tp. 


2tr 


f dY  f d|Ypv  (Y)  exp  [n(r(0)Y*  - |Y|  2/2)/hvW] 

J0  J0  1 


say  H, 


n 


CV.4) 


< 

say  H, 


0 


where  Y exp(j<£>)  is  the  polar  form  of  the  complex  number  Y,  p^  (Y)  is  the 


probability  density  for  the  complex  random  variable  y_,  and  n is  chosen  to 
obtain  a desired  Pp  level.  The  results  of  Sections  II,  III  imply  that 


Y P^QO  = Py(y)/2iT  for  0 < <p  < 2ir 


(V.5) 


where  p (Y)  is  the  probability  density  for  y = |yj . Thus,  the  LRT  (V.4) 


can  be  reduced  to 


/ dYpy(Y)  IoCn|r(0)|Y/hvoW)exp(-nY2/2hvoW) 
0 

say  H, 


> 

— % 

n. 


(V.6) 


< 

say  Hf 
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where  I is  the  zeroth-order  modified  Bessel  function  [3b].  Now.  because 
o 

I is  a monotonically  increasing  function  of  its  argument,  the  LRT  (V.6)  is 
equivalent  to  the  threshold  test 

say 

> 

|r(0)|2  y CV.7) 


say  HQ 

where  y is  chosen  to  achieve  the  desired  false-alarm  probability.  Because 
W = tp  ■*"  has  been  assumed,  (V.7)  may  be  realized  by  matched-filter  envelope 
detection  followed  by  threshold  comparison.  Note  that  the  optimality  of  the 
threshold  test  (V.7)  applies  regardless  of  the  target  character  (specular, 
diffuse  or  semi -rough)  and  whether  or  not  there  is  turbulence  present  along 
the  propagation  path.  The  performance  of  this  test,  viz,  the  ROC,  is  not 
similarly  insensitive  to  the  details  of  the  problem.  We  shall  investigate 
its  behavior  for  specular  and  diffuse  targets. 

V.2  The  Receiver  Operating  Characteristic 

It  turns  out  to  be  relatively  easy  to  obtain  the  receiver  operating 
characteristic  for  the  threshold  test  (V.7)  when  the  target  is  either 
specular  or  diffuse.  Regardless  of  target  and  turbulence  conditions,  the 
IF  complex  envelope  r(t)  is  pure  receiver  noise  under  hypothesis  HQ.  Thus 
for  any  target/turbulence  combination  we  find  that  the  false-alarm  probability 
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IS 


Pp  = Pr  [|n(0)|2  y]  = exp(-ny/2hv0W)  (V.8) 

Equation  (V.8)  allows  us  to  solve  for  the  threshold  level  which  gives  false- 
alarm  rate  Pp,  viz,  the  test 

say  Hp 

> . 

|r(0)|2  -(2hv0W/n)ln  PF  (V.9) 

< 

say  Hfl 

has  the  maximum  P^-value  attainable  with  false-alarm  probability  Pp.  The 
detection  probability  is  given  by 

PD  = Pr  [ |jr  + n(0)  | 2 > - (2hvQW/n)ln  Pp] ; (V.10) 

more  explicit  results  are  developed  below. 

Case  1 Specular  Target 

For  a specular  target  in  the  single-glint  approximation  we  have 
(cf.  (IV. 17)) 

y k P T1/2/  dp'  it2(p’)  Tg(p')exp(2x(pg',  6)  + 2j<),(pg',  O)+j0) 

(V.ll) 

Standard  results  from  Gaussian  detection  theory  [36]  thus  yield  the 
conditional  detection  probability 
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Pr  [|Z  + n(0) | 2 > -(2hv0W/n)ln  Pp|  |x|=YJ 

= (^((nY2/hv0W)1/2,  (-2  In  PF)1/2): 


(V .12} 


where 


OO 

Q^(a,B)  S/  du  u exp  [-(u2  + ot2)/3  IQ(au) 

is  Marcum's  Q-function  [36]- [38],  It  follows  from  (V.ll)  that  |yj  is  a 
lognormal  random  variable.  Indeed,  from  Section  IV,  we  have  that 

(nll|2/hv0W)1/2  = (2CNR°g)1/2  exp(2x(pg' , 6)).  (V.13) 


Because  x(p  0)  is  a Gaussian  random  variable  with  mean  -a  and  variance 
§ X 

2 

a , the  detection  probability  is  therefore 

X 


00 

PD  =/  dx  Px(x)  (^(C2CNR0g)1/2e2x,(-2£n  P{ 


■)L/n  CV.14) 


where  p^Cy)  = (2-rro2)  1 ^ exp[-(x  + o2)2/2o2]. 

Despite  the  plethora  of  communication-theoretic  calculations  that  are 
in  the  literature  of  optical  propagation  through  turbulence,  Eq . (V.14J  has 
yet  to  be  evaluated.  Results  for  Marcum's  Q-function  are  available  [36] -[38], 
so  that  numerical  values  of  the  free-space  detection  probability  (obtained 
by  setting  o2  = 0 in  (V.14)] 

PD°  = Qm((2CNR°  )1/2,  (-2  In  PF)1/2).  (V .15) 


can  be  obtained  without  computation.  In  Fig.  17  we  have  sketched  Pn°,  given 


P°  ( percent) 


I 


by  CV.15),  vs.  CNR°  for  a variety  of  interesting  false  alarm  probabilities. 
From  this  figure  we  see  that  the  example  in  Table  2,  Fig.  11  will  give 
excellent  target  detection  performance,  i.e.,  Pp°  0.98  for  Pp  = 10  ^ in 
a 2dB/km  absorption  atmosphere  out  to  more  than  6 km  path  length. 

Numerical  work  on  (V.14)  is  presently  underway  to  quantify  the  effect 
of  scintillation  on  target  detection  performance.  We  expect,  based  on  error 
probability  calculations  for  conmunication  through  turbulence  [10],  [39], 
that  scintillation  will  require  a significant  increase  in  carrier -to -noise 
ratio  to  achieve  the  desirable  high-PpAdw  -Pp  values. 

Case  2 Diffuse  target 

For  a diffuse  target 

£=  P Tl/1f  dp'  lt2(p’)  Tgfp’)  exp(2x(p’,  6)  + 2j<j)(p' , 0)) 

will  be  a zero-mean  circulo-complex  Gaussian  random  variable  conditioned 
on  knowledge  of  x(p'»  0)  + j d> Cp ' > 0)-  Gaussian  detection  theory  now  gives 
us  [36]  the  conditional  detection  probability 


Pr[|x  + n(0)  | 2 > (-2hv0W/n)  In  Pp  |x(p',0)  + j<j>(p',0)] 

=p  (1  + nPT>2/  dp’  I^Cp')!2  Ts(p')  exp(4x(p,,0))/2hvoW)'1 
F 

$ 

CV.16) 


i 


i 
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Because  sums  of  real -valued  lognormal  random  variables  are  themselves 
approximately  lognormal  random  variables  [40],  [41],  the  random  variable 
u defined  by 

e2u  = nP T*Zy  dp’|it2(p’)|2  Ts(p')e4x(p'^>2hv0WCNRs  (V.17) 

2 2 

may  be  taken  as  Gaussian  with  mean  -a  and  variance  a where  (assuming  Ts 
is  constant  over  the  illumination  region) 

4a2  /dpl'/dp2’  l^t2(Pi’))2  llt2(P2’)|2  expdbCxxCp/-^’,^) 

e = 

c / dp'  llt2(p')|2) 


CV.18) 


4 4 

Equation  (V.18)  exhibits  the  aperture -averaging  effect,  viz,  e4°  = 1 + 
2 

C(e^ax  - 1)  where  £ 'v  1 for  XL/d  < (AL)^2  and  £ ^ d2/XL  for  XL/d  >> 


(XL) 17  L (recall  (IV. 36)). 

Using  the  above  results  we  find  that  the  detection  probability  is 


PD  = f dUpu(U)  PF  (1  + ^s  exp(2U))  1 (V.19) 


_ 2 -1/2  2 2 2 

where  PU(U)  = (2no  ) exp[-(U  + a ) /2a  ].  Equation  (V.19)  has  not  been 

evaluted  in  the  literature;  we  have  it  under  study  at  present.  In  the  absence 

2 

of  turbulence  we  have  a = 0 and  (V.19)  collapses  to  the  well-known  free 
space  result 


p o = p (1  + CNR  °)_1 
rD  *F  s 


(V. 20) 


63 


shown  in  Fig.  18.  This  figure  shows  that  the  example  of  Table  2,  Fig.  12 

-14 

will  give  good  target  detection  performance,  i.e.,  Pg°  _>  0.95  for  Pp  = 10 
in  a 2dB/km  absorption  atmosphere  out  to  more  than  4 km  path  length. 

Comparison  of  Figs.  17  and  18  reveals  that  it  is  significantly  easier 
to  achieve  high-Pp)°/low-Pp  values  for  a glint  target  than  for  a speckle  target; 
for  Pp°  = 0.98,  Pp  = 10  ^ the  glint  target  requires  a 17  dB  carrier-to-noise 
ratio  whereas  the  speckle  target  requires  a 32  dB  carrier-to-noise  ratio. 
Paradoxically,  it  would  seem,  at  low  carrier-to-noise  ratio  it  is  easier  to 
detect  a speckle  target  than  a glint  target;  for  Pp  = 10  and  10  dB  carrier- 
to-noise  ratio  we  find  Pp°  = 0.014  for  a glint  target  and  PD°  = 0.123  for  a 
speckle  target.  The  explanation  of  the  high- CNR  disadvantage  and  the  low-CNR 
advantage  of  the  speckle  target  vs.  the  glint  target  is  as  follows.  When 
CNR  ° = CNR  ° the  average  strength  of  the  randomly  distributed  speckle-target 

& s 

return  equals  the  strength  of  the  non-fluctuating  glint -target  return.  At' 
high  carrier-to-noise  ratios,  speckle-target  detection  probability  is  limited 
primarily  by  the  chance  that  a deep  fade  will  occur  (i.e.,  the  target  return 
is  much  weaker  than  its  mean  value  and  thus  buried  in  the  receiver  noise) , 
hence  the  disadvantage  relative  to  a glint  target.  At  low  carrier-to-noise 
ratios,  speckle -target  detection  probability  is  enhanced  by  the  chance  that 
the  target  return  will  be  much  stronger  than  its  mean  value  and  thus  stand 
out  from  the  receiver  noise;  this  explains  the  performance  advantage  relative 
to  a glint  target. 

VI ; SUMMARY 

In  this  section  we  shall  summarize  the  key  results  of  our  imaging  and 
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-nation  dif^e-target^e^obabi; 


I 


target  detection  analyses.  We  shall  draw  freely  from  Sections  II -V  without 
explicitly  stating  the  conditions  required  for  the  validity  of  the  results, 
and  we  shall  presume  familiarity  with  the  notation  developed  in  those  sections 
(see  Table  3).  We  begin  with  a qualitative  description  of  our  theory. 


Qualitative  Radar  Behavior 

In  all  weather  conditions  for  which  atmospheric  backscatter  can  be 
ignored,  the  radar  receiver's  IF  signal  consists  of  a target  return  plus 
receiver  noise.  The  target  return  may  be  influenced  by  atmospheric  beam 
spreading  (which  can  reduce  the  on-target  irradiance)  and  atmospheric 
scintillation,  as  well  as  the  reflection  properties  of  the  target  itself. 
Because  we  have  considered  a heterodyne  reception  radar,  atmospheric  coherence- 
loss  can  also  affect  the  target-return  IF  signal.  Furthermore,  atmospheric 
absorption  will  introduce  an  overall  exponential  extinction  of  the  target 
return. 

Under  typical  clear-weather  atmospheric  conditions  we  shall  have  d < pQ, 
so  that  atmospheric  beam  spreading  and  receiver  coherence-loss  may  be  neg- 
lected. Scintillation,  however,  will  play  a significant  role  in  both  imaging 
and  target-detection  applications.  In  the  former  instance,  scintillation  will 
limit  achievable  image  signal -to-noise  ratio  for  a glint  target  observed  under 
high  carrier-to-noise  ratio  conditions;  for  a speckle  target,  scintillation 
may  reduce  image  signal -to-noise  ratio  significantly  below  the  unity  level 
set  by  the  target  fluctuations.  Multiframe  averaging  will  be  needed  to  over- 
come the  scintillation  limit  on  image  signal -to -noise  ratio. 


5" 
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TABLE  3 

SUMMARY  OF  NOTATION 


Radar  Parameters 


P.j,  = transmitter  peak  power 

d = optics  diameter 

A = 10.6im  wavelength 

hvQ  = photon  energy  at  10.6pn  wavelength 

n = detector  quantum  efficiency 

2W  = unilateral  IF  bandwidth 

= normalized  target -plane  illuminator  beam 
Atmospheric  Parameters 


X 

<P 


a 


turbulence- induced  log-amplitude  fluctuation 

turbulence- induced  phase  fluctuation 

field  coherence  length 

path  length 

log-amplitude  variance 

absorption  coefficient 

aperture -averaging  factor 

aperture -averaged  log -amplitude  variance 


Target  Parameters 


T = diffuse-reflection  coefficient 

" 5 
T 


specular-reflection  coefficient 
Aj,  = target  area 

p* (Ajf^jfy)  = target  bidirectional  reflectance 
T = average  diffuse-reflection  intensity  coefficient 
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! In  single-pulse  target  detection,  target  speckle  and/or  scintillation 

dictate  that  receiver  carrier -to -noise  ratio  must  be  increased  (over  that 
i needed  for  a free-space  propagation/glint -target  system)  to  achieve  the 

desirable  combination  of  a high  detection  probability  at  a low  false-alarm 
rate. 

IF  Signal  Complex  Envelope 

j Assuming  d < pQ  and  a target  reflection  coefficient  with  diffuse  and 

single-glint  components,  we  have  that  the  IF  signal  complex  envelope  satisfies 

r(t)  = £ + n(t)  (VI.  1) 

where 

Z^PT 1/2 f dp'  h2(~p,)  exp(2x(p\0)  + 2j <t>(p'  ,0)  -aL) 

♦ P T1/2/  dp'  it2(p')  Tgfp' ) exp(2x(pg,,0)  + 2jd>(pg’  ,6)  + je-oL), 

(VI. 2) 

is  the  target  return,  and  n(t) , a zero-mean  circulo-complex  Gaussian  noise 
process  with  spectral  density 

Snn(f)  = hvo/n  (VI.  3) 

i 

is  the  heterodyne -receiver  shot  noise. 

Carrier-to-Noise  Ratio 

The  carrier- to -noise  ratio  for  the  above  IF  signal  envelope, 
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CNR  = <|z|2>  / <|n|2>,  (VI .4) 

is  given  by 

CNR  = CNR  + CNR  (VI. 5) 

8 s 

in  terms  of  its  glint  contribution  (Fig.  11) 

2 

CNRg  = (nPT/2hvoW)|^.2CXLfT)|2  A2  Aj.  p’ (A;fT;-fT)e4aX  "2aL'  (VI. 6) 
and  its  speckle  contribution  (Fig.  12) 

2 

CNRS  = (nPT/2hvQW)  A2  Ts(ALfT) f dp’  |£t2(p’)|2  e4°X  '2aL  (VI. 7) 


Image  Signal-to-Noise  Ratio 

The  image  signal-to-noise  ratio, 


(< | r (0) |2>  - 2hv  W/n)2 

SNR  = =2 (VI.  8) 

Var(|r(0)|2) 


satisfies 


SNR  = 2^5/2 , (VI.  9) 

1 + CNR/2SNRSAT  + (2CNR)_i 

where 

SNRjjat  h <|Z|2>2  / Var(|Z|2)  (VI.  10) 

is  the  saturation  signal-to-noise  ratio  (due  to  target  and  turbulence 
fluctuations).  For  CNR  >5,  (VI. 9)  reduces  to  the  universal  curve 


1 

i 

f 

I 

I 


i 
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SNR  y CNR/2SNRSAT 

'X,  - ■ y 

SNRsat  1 + CNR/2SNRsat 


which  shows  (Fig.  9)  that  increasing  CNR  beyond  2SNRgAT  does  little  to  improve 
image  signal-to-noise  ratio. 

The  saturation  signal-to-noise  ratio,  (VI. 10),  has  the  following 
behavior: 

i)  No  turbulence,  specular  target 


SNR' 


gSAT 


ii)  No  turbulence,  diffuse  target 


SNR°sSAT  = 1 


iii)  No  turbulence,  semi -rough  target 


SNR°SAT  = 1 + CNR0g/GNR0s 
iv)  Turbulence,  specular  target  (Fig.  13) 

SNRgSAT  X _ 

v)  Turbulence,  diffuse  target  (Fig.  15) 

2 

SNRsSAT  = [1  + 2te16ax  - D5]-1 


(VI. 12) 


(VI. 13) 


(VI. 14) 


(VI. 15) 


(VI. 16) 


Likelihood -Ratio  Test 

The  decision  rule  which  maximizes  the  conditional  probability  you  say 
a target  is  present  given  there  is  one  present,  subject  to  the  constraint 
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that  the  conditional  probability  you  say  a target  is  present  given  there  is 
no  target  there,  is  the  threshold  test 


say 

> 

|r(0)|2  -(2hv0W/n)to  PF 

< 

say  HQ 

Receiver  Operating  Characteristic 

The  threshold  test  (VI. 17)  has  false  alarm  probability  Pp 
probability  given  as  follows: 

i)  No  turbulence,  specular  target  (Fig.  17) 

PD°  = (^((2CNR°g)1/2,  (-2ta  PF)1/2) 

where  is  Marcum's  Q function 

ii)  No  turbulence,  diffuse  target  (Fig.  18) 

p o _ p (1  + CNR0!"1 
D *F 

iii)  Turbulence,  specular  target 

00 

PD  = / dx  Px(x)  <^((2CNR0g)1/2  e2X,(-2£n  Pp)1/2) 
-00 

where  Py(x)  = (2ito2)'1/2  expf-(x  + a2)  /2a2] 


(VI. 17) 


and  detection 


(VI. 18) 


(VI. 19) 


(VI. 20) 
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